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We formulate Dirac fermions on a (l+l)-dimensional lat- 
tice based on a Hamiltonian formalism. The species doubling 
problem of the lattice fermion is resolved by introducing hop- 
ping interactions that mix left- and right-handed fermions 
around the momentum boundary. Approximate chiral sym- 
metry is realized on the lattice. The deviation of the fermion 
propagator from the continuum one is small. 



In contrast with the great success of lattice gauge the- 
ory, lattice fermions remain a long-standing problem. 
Naive discretization causes the well-known species dou- 
bling problem [1]. The problem originates in the fact 
that the kinetic term of the fermion is proportional to 
the first-order derivative in real space. This means that 
the Fourier transform of the kinetic term is proportional 
to the momentum. Since an odd function cannot be peri- 
odic and have one zero, the fermion propagator is forced 
to have additional pole(s) on the lattice. The situation 
does not change regardless of how the lattice spacing is 
reduced as long as the space-time derivative is modeled 
as a naive difference. 

Many attempts have been made to fix the doubling 
problem [1-7]. Wilson removed doublers at low energy 
by introducing an interaction that mixes left- and right- 
handed fermions [1]. The interaction is not a result of 
naive discretization and therefore has no counterpart in 
continuum theory. However, unwanted degeneracy per- 
sist at high energy and chiral symmetry is explicitly bro- 
ken. To fix these problems, Kaplan modified Wilson's 
fermion by introducing an extra dimension [4] . Kaplan's 
fermion has an approximate chiral symmetry if the lattice 
size of the extra dimension is large. The fermion is useful 
for calculating physical quantities related to dynamical 
breaking of chiral symmetry. (See Ref. [8], for exam- 
ple.) However, the cost of numerical calculations based 
on Kaplan's fermion is not cheap. If we find a method to 
perform such calculations without the extra dimension, 
calculation time decreases largely and a deeper under- 
standing of quantum field theory becomes possible. 

In addition to the doubling problem, the lattice 
fermion has another serious problem. The fermion prop- 
agator defined on a lattice deviates from the continuum 
one even if the doublers are removed with the existing 
techniques such as Kaplan's fermion [5]. As a result, it 
may cause errors in numerical calculations. To remove 
this uncertainty, we need to modify the discretized prop- 
agator somehow so that it is close to the continuum one 



as far as possible. Such deviation of propagators becomes 
critical especially when supersymmetry is considered on 
a lattice because deviation of fermion and boson propa- 
gators gives wrong values for loop integrals (sums). For 
example, the zero-point energy does not cancel between 
fermions and bosons if the propagators deviate from the 
continuum ones. Accurate discretization of the propa- 
gators is a necessary condition for maintaining lattice 
supersymmetry. 

In general, the extra dimension can be expressed as 
hopping interactions in a lower-dimensional system. Ka- 
plan's fermion is a formulation with an extra dimension, 
so there must be a corresponding Hamiltonian with no 
extra dimension. Also, the shape of the fermion propa- 
gator can be improved with hopping interactions. The 
Runge-Kutta method for differential equations is an ex- 
ample of such an improvement. 

In this paper, based on a Hamiltonian formalism, we 
introduce ultralocal hopping interactions to remove dou- 
blers and improve momentum dependence of fermion en- 
ergy. (The word "ultralocal" means that fermion hopping 
is restricted to a finite range on a real-space lattice [9].) 
From knowledge of the continuum theory, we know the 
correct momentum dependence of the energy. We start 
from momentum space and go back to real space by way 
of discrete Fourier transform. A real-space Hamiltonian 
is necessary to construct gauge theory. The method is a 
hybrid of the Wilson [1] and SLAC [3] approaches. 

First, let us consider a free Dirac fermion on a (1+1)- 
dimensional Hamiltonian lattice. Time is continuous and 
space is discrete. The length of the spatial lattice is Na, 
where a is a lattice spacing and the number of sites N is 
assumed to be even. According to the continuum theory, 
the Dirac fermion should be described by the following 
Hamiltonian in momentum space 



H = 



1 



N/2 



(1) 



l=-N/2+l 



where I is an index for momentum. Q and Q are discrete 
Fourier transform of real-space two-component fermion 
operators 
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0, (2) 



l=-N/2+l 



and V'n = "0^7" for rt = 1,2, . . . , N, respectively. ^„ and 
rjn satisfy 
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and other anticommutators are zero. The gamma matri- 
ces are 
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Periodic boundary conditions are assumed in real space, 

2-kI 



Pi = 



N 



We try to create a real-space Hamiltonian with no dou- 

blcr that reproduces Eq. (1). To find rcal-spacc repre- 
sentation of pi , let us consider the following function 



M 

E 



2(-l)° 
a 



■ sm ap. 
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In the limit M ^ oo, the function (3) goes to p. The 
function (3) necessarily has a node at the boundary 
p = ±7r because it has a periodicity of 27r (see Fig. 1). 
The node is the cause of the doubling problem. The dou- 
bler remains as a singularity at the boundary even if the 
limit M — » 00 is taken (by "doubler" we mean unwanted 
energy degeneracy that is not contained in the continuum 
theory). Anyway, the parameter M needs to be small 
for practical formulation because M corresponds to the 
maximum distance of fermion hopping in real space. The 
limit M oo needs the infinite lattice. 
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In addition to the doubler modes around the boundary, 
the function (3) has another degeneracy. The function os- 
cillates around p and has local minima if the summation 
is truncated with a small M. In Fourier analysis, it is 
called the Gibbs phenomenon, which occurs if a function 
to be expanded has a singularity [10] . The oscillation can 
be removed by replacing Eq. (3) with the function 



M 



ap, 



(4) 
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where 



Sa = Fr 
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Fa is called the Lanczos factor [10]. As shown in 
Fig. 1, the factor almost removes the oscillation of Eq. 
(3). However, the doubler modes around the boundary 
still remain. We are going to remove them by a trick 
with hopping interactions. Let us consider the following 
momentum-space Hamiltonian: 

N/2 

H= Yl i^iCn'Ci + mCiCi), (5) 

l=-N/2+l 

where s; = s{pi)/a and m is fermion mass. We write the 
Hamiltonian in the matrix form: 



N/2 



l=-N/2+l 



(6) 



We introduce an interaction cj that mixes left- and right- 
handed fermions 



H 



N/2 

l=-N/2+l ^ 



Si 



m + ci 



m + ci -si 



(7) 



where c; are assumed to be nonzero only for jZj ~ N/2. As 
shown later, si and c/ are expressed as ultralocal hopping 
interactions in real space. The Hamiltonian (7) can be 
diagonalized for each I. 
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FIG. 1. The solid line plots the correct energy p 
of one-particle states from the continuum theory. The 

dot-dashed line plots the function (3) for M = 5. The dashed 
line plots the function s(p) modified with the Lanczos fac- 
tor for M = 5. The oscillation of Eq. (3) has almost been 
removed. 



where ki = sj + {m + ci)^ arc energies of one particle 
states and C/ are transformed variables. Although bet- 
ter solutions may be found than c; shown here, we give 
precedence to simplicity over accuracy in this paper. We 
understand that our strategy is successful if the proper- 
tics of the continuum Dirac fermion are approximately 
reproduced. 

Consider the function 



2 



—u{p — tt)^ + = -^^ + lim Ca cos(ap), (9) 

a=l 

where u and v are some positive real numbers. The 
Fourier coefficients Ca are 

Ca = Fai-ir^ [-/cos (a/) + i sin (a.[^) 

The function (9) has a peak at p = ±7r and zeros at 
p = ±7r =F As before, the infinite M cannot be 

reaUzed on a finite lattice. For a finite M, we define the 
function 

M 

c(p) = ^ + VC«cos(ap). (10) 

When M is finite, c(p) does not reproduce the function 
(9) correctly. However, it does not matter in this consid- 
eration. We just want to use c(p) to remove the doubler 
modes. We do not need to care about the original shape 
of the function (9). The parameters u and v are adjusted 
so that the unnecessary degenerate modes around the 
boundary become non-degenerate normal modes. Fig- 
ure 2 shows how the doubler modes of s(p) are removed 
using the function c(p). The function s^(p) has two 
hemlines near the boundary because it has to vanish at 
■p = ±7r. Wc want to raise the both ends by adding a 
packet function {p) to {p) . Let us consider a function 
k{p) = \/ s'^{p) -\- c^ip). If we choose M = 5, u = 130, 
and V = 8.4, the function k'^ip) agrees well with in the 
fundamental region < tt except for a small deviation 
around momentum \p\ ~ 2.3. 
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FIG. 2. The solid line plots the correct energy squared p^. 
The dashed line plots the function s^{p) for M = 5 with the 
Lanczos factor. The dotted line plots the function c'^(p) for 
M = 5, u = 130, and v = 8.4 with the Lanczos factor. The 
dot-dashed line plots the function k^{p). The doubler modes 
around the boundary |p| = 7r has been removed with c^(p). 

In Fig. 3, the functions ±p, ±s{p), and ±k{p) are com- 
pared. The hmction ±k{p) corresponds to one-particle 
states given by the Hamiltonian (7) with m = 0. The 
function ±fc(p) agrees well with the correct energy ±p 
from the continuum theory in the fundamental region 
\p\ < TT. The Hamiltonian (7) approximately reproduces 
the continuum theory for the (l+l)-dimensional Dirac 
spinor without doubler modes if we identify Ci = c{pi)/a 
and ki — k{pi)/a. 
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FIG. 3. The solid line plots the correct energy ±p from the 
continuum theory. The dashed line plots the function ±s(p) 
for M = 5 with the Lanczos factor. The dot-dashed line plots 
the function ±fc(p) for M — 5, u = 130, and v = 8.4 with the 
Lanczos factor. The function ±fc(p) almost agrees with ±p in 
the fundamental region |p| < n except for a small deviation 
around momentum |p| ~ 2.3. 

In the new basis that diagonalizes Eq. (7) with m = 0, 
75 is transformed into 

, _ si + ki f si -cA , . 

kf + Slkl\-Cl -SI J 

for i > and 

, _ -si + ki (si ci \ , . 

for / < 0. For I = 0, the Hamiltonian (7) is diagonal with 
degenerate zero energy, so 75 is not transformed. 
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Figure 4 shows the diagonal and ofF-diagonal matrix 
elements of the transformed 75. The diagonal (1,1) ele- 
ment of Eq. (11) is almost unity at low and intermediate 
energy pi < 2.3 and deviates from unity at pi > 2.3. The 
off-diagonal (1,2) element of Eq. (11) oscillates around 
zero at pi < 2.3 and becomes unity at pi = w. At low en- 
ergy, the deviation of the off-diagonal elements from zero 
is not large and becomes smaller as the parameter M 
increases. At pi < 2.3, the transformed left- and right- 
handed fcrmions have approximately the correct chiral 
charges 1 and —1, respectively. The low-energy Hamil- 
tonian (7) has approximate chiral symmetry because the 
commutation relation between the Hamiltonian (7) and 
chiral charge defined with 75 is almost zero for small I. 
The errors associated with chiral symmetry can be im- 
proved in a systematic way by increasing M. The value 
used here for the parameter M is sufficiently small and 
does not deny application of the model to actual numer- 
ical analysis with a computer. 
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fermions in Eq. (13). The contimmm limit a ^ is 
taken with the parameter M fixed. The method can be 
extended to higher-dimensional cases by taking care of 
the doubler of each direction in the same way. 

For Monte Carlo analysis, it is useful to rewrite our 
Hamiltonian into a Euclidean formulation. The Eu- 
clidean action for the Dirac spinor is given by 



M 2 



n ^ c«=l (U=l 



+ (m+^ )^„^„!>, (14) 



where n indicates a site on a two dimensional Euclidean 
lattice and is a unit one-site vector in the jj, direc- 
tion. The Euclidean gamma matrices satisfy jj^ = 7^, 
and {7p,7i/} = 2S^^ and -tp = The definition of a 

lattice gauge theory based on the action (14) is ambigu- 
ous. An exponentiated gauge field can take any path 
between two ends. The most natural choice for gauge 
field is the shortest path that links the two ends because 
fermion hopping is parallc;! to one of the directions /i. 

In this paper, we have constructed doubler-free Hamil- 
tonian and Euclidean action for Dirac fermions on a 
(l+l)-dimcnsional lattice. To realize approximate chi- 
ral symmetry at low energy, explicit breaking of chiral 
symmetry has been compressed to around the momen- 
tum boundary. In future works, it should be precisely 
checked if insertion of gauge interactions affects chiral 
properties. 

I am indebted to Michael Creutz for valuable discus- 
sions. I also thank Jun-ichi Noaki for informative discus- 
sions. This research was supported in part by RIKEN. 
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FIG. 4. The matrix elements of Eq. (11) are plotted as 
functions of the momentum pi for a = 1, M — 5, u — 130, 
and V = 8.4 with the Lanczos factor. For a lattice size N = 40, 
the functions are plotted in the right-half plane of momentum 
space pi > 0. The circles are plots of the diagonal (1,1) el- 
ement of Eq. (11). The crosses are plots of minus of the 
off-diagonal (1,2) clement of Eq. (11). The solid and dashed 
lines are plotted to guide the eyes. 

To obtain the real-space Hamiltonian for the Dirac 
spinor, we substitute the discrete Fourier transform (2) 
into Eq. (7) 

n=l ^ a=l 



t 



+ |m+^)VnV'n^ (13) 



As usual, gauge symmetry can be implemented by in- 
serting exponentiated gauge fields between the hopping 
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